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Various sum rules for multiple giant dipole resonance states are derived. For the triple 
»— p' giant dipole resonance states, the energy-weighted sum of the transition strengths requires a 

model to be related to those of the single and double giant dipole resonance states. It is also 
^ , shown that the non-diagonal matrix elements of the double commutator between the dipole 

operator and the nuclear Hamiltonian give useful identities for the excitation energy and 
transition strength of each excited state. Using those identities, the relationship between 
width of the single dipole state and those of the multiple ones is qualitatively discussed. 
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It is well known that sum rules play an important role in a wide range of 
physics. ' On the one hand, if an experiment shows breaking of the appropriate 
sum rule, the basic framework of the model should be improved. For example, 
the enhancement of the Thomas- Reiche-Kuhn(TRK) sum rule value requires meson- 
exchange currents. 1 )' 2 ) The quenching of the Coulomb 3 ) and Ikeda 4 ) sum rule values 
indicates that nucleon-degrees of freedom are not enough for understanding of the 
relevant nuclear phenomena. On the other hand, when we make some approxima- 
tions in calculations, we should take care of the fundamental sum rules. From this 
point of view, the sum rules in the random phase approximation have been explored 
by many authors in non-relativistic 5 ) and relativistic 6 ) frameworks. 

In nuclear physics, a typical example of the sum rules is obtained by the rela- 
tionship 



X 



[D,[H,D]] = C, (1) 

where H,D, and C denote the Hamiltonian of the system, the excitation operator, 
and a c-number, respectively. When we take the ground-state expectation value of 
the both sides, we have the model-independent sum rule, 

2j> n |(n|D| )\ 2 = C, (2) 



where u n represents the excitation energy of the eigen state of the Hamiltonian |n) 
from the ground state | ). For example, the TRK sum rule for dipole excitations, 



(3) 
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is obtained for 

i=l i=l 

assuming its double commutator with the nuclear interactions to be zero. Recently 
in the same approximation as for the TRK sum rule, another sum rule has been 
derived, 7 ) 

^u; n \(n\D 2 \ )| 2 = 45 (l)5i(l), (5) 

n 

where Sx(l) stands for the TRK sum rule value and Sq(1) is the non-energy weighted 
sum of the dipole transition strengths, 

S (l) = J2\(n\D\ )| 2 . (6) 

n 

The new sum rule Eq.© is related to the double giant dipole resonance(DGDR) 
states which have been recently observed through relativistic heavy ion reactions. 8 ) 
It implies that the energy-weighted sum of the transition strengths for the DGDR is 
determined by the excitation energies and strengths of the single giant dipole reso- 
nance(SGDR) states. This relationship is useful both for the analysis of experimental 
data and for making models of the DGDR, since the SGDR is well known. 

The purpose of the present paper is to derive other sum rules for multiple dipole 
states on the basis of the same assumption as for the TRK sum rule. In particular, 
we will discuss the sum rules for the operator D 3 , which is related to the triple giant 
dipole resonance(TGDR) states. We will also show that the non-diagonal matrix 
elements of Eq.(^Q) provide us with useful identities to understand the structure of 
the multiple giant resonance states. They give a constraint not only on the sum of 
the matrix elements, but also on each energy- weighted matrix element. 

First we discuss sum rules for the operator D 3 . For convenience, let us define 
the notation for the sum of the matrix elements, 

S k {m) = Y,u k mn \(mn\D m \ )| 2 , (7) 

n 

where oo mn denotes the excitation energy of the ra-th state of the m-times SGDR, 
\mn). In this notation, the TRK sum rule Eq.@ and the new sum rule Eq.(J5J) 
mentioned above are expressed as 

5i(l) = J>i„|(ln|D| )| 2 , (8) 

n 

S 1 (2) = Y,"2n\(2n\D 2 \ )| 2 . (9) 

n 

The above two sum rules are calculated with the use of the closure property of 
the intermediate states. For example, 51(2) is expressed as 



Si(2) = - V( \[D 2 ,H}\2n)(2n\D 2 \ ) + ( | D 2 | 2n ){ 2n | [H, D 2 ] | ) . (10) 
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This equation is rewritten with the closure property as 

Si(2) = l< \[D 2 ,[H,D 2 }]\ ). (11) 

By calculating explicitly the double commutator, we obtain Eq.Q. 

Unfortunately, for Si (3) we can not use the same procedure, because the final 
states excited by D 3 are not only |3n), but also \ln). In the case of Ea. (|lU|) . D 2 
has also a matrix element between the ground states, but it does not contribute 
to the energy-weighted sum and therefore, we could use the closure property. In 
other words, the sum rule value obtained by the double commutator of with 
Hamiltonian is not the value of Si (3). Keeping this fact in mind, we calculate 

S 1 (f) = Y,"f\(f\D 3 \ )| 2 

/ 

= \{ \ [D 3 ,[H,D 3 ]] \ ) = 9S 1 (1)(S 2 (1) + S (2)), (12) 

where we have used 

( \D*\ ) = S 2 (1)+S (2). (13) 

Thus the value of S\(f) is fixed by the TRK sum rule values and the transition 
strengths of the SGDR and DGDR. In the above equation, however, |/) should be 
either |ln) or |3n). Fig. 1 shows the transitions to be included in Si(/). In addition 
to (a) corresponding to Si (3), four kinds of transitions contribute to Si(/). 




o 



(a) (b) 



n i 



(c) 



(d) 



(e) 



Fig. 1. Transitions induced by D 



Among the transitions in Fig. 1, the sums of the strengths in (c), (d), and (e) 
are calculated individually. The contribution (c) is written as 

(c) 



i',n",n m 



./( \D\\n){\n\D\ ){ \ D \ ln')( In' | D \ )( | D \ In" )(n" \ D\ ) 



Si(l)S 2 (l) 



(14) 



4 



H. Kurasawa, T. Suzuki and P. F. Bortignon 



The sums of (d) and (e), on the other hand, are described as 

(d) = (e) = u ln »( \D\ln){ln\D\2n'){2n'\D\ln") 

n,n',n",n"' 

x (ln"\D\ )( \D\ln"')(ln'"\D\ ) 

= \{ \ [D\[H,D}\ \ )5 (1)-5 1 (1)5 2 (1)=2S 2 (1)5 1 (1). (15) 

From Eas.(|12|). Q14|) and (|15|) . we have the sum of (a) and (b) 

Si (3) + ^2 ui n »{ | D I ln)( In | D 1 2n'){2n' \ D \ In" ) 

x ( In" | D | 2n w ) ( 2n'" \ D \ In"" ) ( In"" \D\ ) 
= 4S 1 (1)S 2 (1)+9S 1 (1)S (2). (16) 

Thus, although several sum rules are obtained for D 3 , we can not obtain the re- 
lationship of Si (3) to the energy- weighted and non-energy- weighted sum for the 
strengths of the SGDR and DGDR. The second term of the left hand side in the 
above equation, coming from (b), can not be expressed in terms of Sfc(m). 

In order to separate the contribution Si (3) from (b), we need an assumption. 
For example, if we assume the folding model where the DGDR states are composed 
of the dipole bosons |ln), 

\2n) = \ln',ln") (ri >n"), (17) 

then we have 

(b)=4Si(l)S 2 (l). (18) 
Here, we have used the equation of the folding model, 

^2{ln\D\2n"){2n" \D\ln ) = <WS (1) + ( | D | In ){ In | D | ), (19) 

n" 

which is obtained from the fact that 

( In | D | In', In" ) = y/2( \D\ln )S nn >5 n > n » + ( | D \ In" }<W(1 - <5„' n ") 

+ ( |£|ln')ff„„//(l-<W). (20) 

Eqs.(|inj) and (JTHJ) yield 

Si (3) = 9Si(l)S (2). (21) 

We note that Passos et al. 9 ) have defined TGDR states which are excited by the 
operator 

q( 3 ) = D 3 — D - j . (22) 
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This operator is determined with the requirement that the doorway state of the 
TGDR is orthogonal to the one of the SGDR, 

J2( \D\ln)(ln\0 (3) \ ) = 0. (23) 

n 

Their TGDR states defined in this way are different from the present |3n) which 
satisfies (ln\3n') = 0. 

The sum rules discussed so far are related to the diagonal matrix elements of 
Eq.Q. Next we will show that if we calculate its non-diagonal matrix elements, we 
can obtain other useful identities. 

The equation 

(m + l,n | [D, [H,D]] \m - l,ri) = (24) 

provides us with 

+ ^m-i,n>) ^( m + l,n\D\ mn")(mn" \D\m — l,n') 

n" 

= 2 u> mn "(m + 1, n | D \ mn" )(mn" \ D\m — 1, n'). (25) 

n" 

For m=l and 2, for example, we have the relationships, 

uJ 2n {2n\D 2 \ ) =2j2 UJ i n i(2n\D\ln f )(ln'\D\ ), (26) 

n' 

uJ 3n (3n\D 3 \ ) =3^2uJ lnl (3n\D 2 \ln')(ln \ D\ ). (27) 

n' 

In deriving the last equation, we have used Ea.(|26[). together with Ea.l|25|l. According 
to these identities, it is possible to make the following comments. 

First, the excitation energy u>2n is determined by io\ n and transition strengths. 
This fact is also true for other uj mn . The sum rule Eq.(JSJ) is the relationship between 
the sum of the energy-weighted strengths for the DGDR and those of the SGDR, 
while Ea. (|26|) provides the relationship of each energy- weighted strength of \2n) to 
uj\ n and transition strengths, 

uj 2n { \D 2 \2n)(2n\D 2 \ ) =2^w ln /{ | D 2 \ 2n)( 2n \ D \ In )( In \ D \ ). (28) 

n' 

In using the above equation, the energy-weighted sum of the DGDR is expressed 
in the form: 

5i(2) = 2j^wi n /< \D 2 \2n){2n\D\ln'){ln'\D\ ). (29) 

n,n' 

Since the right-hand side of this equation is expressed as 

v.h.s. = 2^2uj ln ( \D 3 \ln)(ln\D\ } - 2S (l)Si(l), (30) 

n 
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we have another sum rule, 

|£> 3 |ln)(ln|£>| ) = 3S (l)Si(l). (31) 

n 

Of course, it is possible to derive the last sum rule by calculating the double commu- 
tator, as was done in Ea.()15|). In this case, the sum rule of the DGDR, Eq.(JSJ), can 
be obtained from Ea. (|26j) . We note also that Eq. (|15|) is easily obtained in employing 
Eq.((2EI>. 

Second, if there is a single collective state \m = 1} with the excitation energy u>i, 
then Eq.(J2ni) yields w% n = 2cji, and hence, owing to Eg . (|23)) . we have uj mn = muj\. 
More approximately, Ea. (|26j) implies that if uj\ n ~ a), we should have u>2n ~ 2tD, 
where ui represents the mean energy of the SGDR. This fact is expressed by rewriting 
Eq.(|inj) as 

uj2n(2n\D 2 \ ) = 2^{ui nl - uj){2n\D\ln'){ln' \D\ ) + 2cu( 2n \ D 2 \ ). (32) 

n' 

Thus, if the width of the SGDR is narrow, the width of the DGDR is also expected 
to be narrow. The conclusion on other m-GDR is also the same. In other word, if the 
width of the SGDR is narrow, interactions satisfying [D, [V, D]] = may not play a 
role to make the width of the m-GDR broader. On the one hand, the diagonal matrix 
element of Eq.Q shows that there should be the same dipole strength on any state. 
This fact is usually called Brink's hypothesis. 8 ^ On the other hand, its off-diagonal 
matrix elements provide a constraint on the distribution of the strengths. 

The constraint on the width of the DGDR due to the SGDR may be seen more 
clearly in the following identity which is obtained from Ea, (|26|) . 

Y,(u2n-2u) 2 \(2n\D 2 \ )| 2 

n 

= 4:^2(u ln -Q)(u ln > -Q){ \D\ln)(ln\D 2 \ln')(ln \D\ ), (33) 

n,n' 

where the mean energy of the SGDR is defined by to = S\(l)/ So(l). 

In the folding model, the above equation yields the relationship between the 
variances of the DGDR and SGDR strength functions. Ea. (|19|) gives 

So(2) = 2,S 2 (l), (34) 

and 

Y,U2n\(2n\D 2 \ )\ 2 /S q (2)=2Cj. (35) 

n 

The straightforward calculation of Ea. (|33|) . together with the above two equations 
and Ea. ()19|) . provides us with the well-known result, 7 )' 8 - 1 



<7 2 = V2(Jl, 



(36) 



Sum Rules of the Multiple Giant Dipole States 



7 



where o~\ and a"2 denote the variances of the SGDR and DGDR, respectively, 

al = Y,(uin-0) 2 \(ln\D\ )\ 2 /S (l), (37) 

n 

al = Y,{^n-2ujf\{2n\D 2 \ )\ 2 /S (2). (38) 

n 

It should be noted that the broadening by the factor v2 in Ea.(|36|) is not due to 
nuclear interactions. This is just a stochastic factor given for the two-dimensional 
probability distribution, when there is no correlation between the two stochastic 
variables. In the folding model, the stochastic variable is the excitation energies of 
the bosons, and the probability function is related to the transition strength function. 
Since the covariance in the two boson sets is zero in the folding model, the only first 
term of the right hand side in Ea.(|19|) contributes to the right hand side of Eq. ()33|) . 
Eas. ()34j) and (|35|) are also derived from a stochastic point of view. In the folding 
model, we have a m = ^fma\. 

In summary, various sum rules for multiple giant dipole states are derived by 
assuming, in the same way as for the Thomas-Reiche-Kuhn sum rule, that the dou- 
ble commutator of the dipole operator with the nuclear Hamiltonian is a constant. 
Nuclear interactions depending only on nucleon coordinates, exchange forces with 
zero-range and spin-orbit forces satisfy this assumption. It has been discussed that 
the sum rule for the triple giant dipole resonance can not be obtained without further 
assumptions. 

The sum rules are derived from the diagonal matrix elements of the double 
commutator. We have shown that non-diagonal matrix elements also provide us 
with useful identities, which yield constraints on the excitation energy and transition 
strength of each nuclear state. By using those identities, we can discuss qualitatively 
the relationship between the width of the single giant dipole resonance and those 
of the multiple ones. If the width of the single dipole giant resonance is narrow, 
then that of the multiple one is expected to be also narrow. For more quantitative 
understanding of excited states, and effects of exchange forces with finite range or 
velocity-dependent forces, of course, more elaborate calculations are required. 

Finally, we note that in fact, throughout this paper, we did not use the explicit 
form of the dipole operator, but assumed only the double commutator of the exci- 
tation operator with the Hamiltonian to be constant. Therefore, the present study 
may be useful not only in nuclear physics, but also in discussions of other quantum 
systems like atomic clusters and Bose-Einstein condensation. Moreover, collective 
excitations other than the dipole one may be discussed in the same way. The new 
identities, derived from the non-diagonal matrix elements of the double commutator, 
may be useful for discussions on anharmonicity of collective motions in many-body 
systems. 

References 

1) T. Suzuki, Ann. de Phys. 9, 535 (1984), and references therein. 

2) A. Bohr, and B. R. Mottelson, Nuclear structure, vol.2 (Benjamin, New York, 1975). 

3) H. Kurasawa and T. Suzuki, Nucl. Phys. A490, 571 (1988); Prog, Theor. Phys. 86, 773 
(1991); BgNS Trans. 7, 3 (2002). 



8 



H. Kurasawa, T. Suzuki and P. F. Bortignon 



4) H. Kurasawa, T. Suzuki, and Nguyen Van Giai, Phys. Rev. Lett. 91, 062501, and references 
therein. 

5) D. J. Thouless, Nucl. Phys. 22, 78 (1961); G. F. Bertsch and S. F. Tsai, Phys. Rep. Phys. 
Lett. 18c, 125 (1975). 

6) H. Kurasawa and T. Suzuki, Phys. Rev. C 69, 014306 (2004). 

7) H. Kurasawa and T. Suzuki, Nucl. Phys. A597, 374 (1996) 

8) T. Aumann, P. F. Bortignon, and H. Emling, Ann. Rev. Nucl. Part. Sci. 48, 351 (1998); 
Ph.Chomaz and N. Frascaria, Phys. Rep. 252, 275 (1995); C.A. Bertulani and V.Yu. Pono- 
marev, Phys. Rep. 321, 139 (1999). 

9) E. J. V. de Passos, M. S. Hussein, L. F. Canto, and B, V. Carlson, Phys. Rev. C 65, 034326 
(2002). 



